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: Abstract 



J> ■ The perturbation theory with a variational basis is constructed and analyzed. 

■ The generalized Gaussian effective potential is introduced and evaluated up 

■ to the second order for selfinteracting scalar fields in one and two spatial 

Q^ ■ dimensions. The problem of the renormalization of the mass is discussed 

in details. Thermal corrections are incorporated. The comparison between 
Oh! 

1> ■ the finite temperature generalized Gaussian effective potential and the finite 

temperature effective potential is critically analyzed. The phenomenon of the 
restoration at high temperature of the symmetry broken at zero temperature 
is discussed. 

PACS numbers: ll.lO.-z, ll.lO.Kk, ll.15.Tk, ll.lO.Wx 



X 



1 



Typeset using REVTgX 



I. INTRODUCTION 



In the recent years the variational Gaussian approximation has played an important role 
in the non perturbative study of quantum field theories. In particular, to investigate the 
spontaneous symmetry breaking phenomenon in scalar quantum field theories it has been 
introduced the Gaussian effective potential |jT|]. The main disadvantage of the variational 
approach is the absence of the control of the approximation. Moreover, in quantum field 
theories the presence of ultraviolet divergences often makes useless the variational calcula- 
tions. 

The aim of this paper is to develop a variational scheme in scalar quantum field theories 
which allows to evaluate in a systematic manner the corrections to the Gaussian approxima- 
tion and, at the same time, to keep under control the ultraviolet divergences. To this end, we 
shall construct a perturbation theory with a variational basis. The method we shall follow 
is widely used in many-body theory where it is known as the method of correlated basis 
functions [Q]. Using the variational basis we construct a vacuum state \Q) which is adiabat- 
ically connected to the Gaussian trial vacuum. Whereupon, we introduce the generalized 
Gaussian effective potential Vg{4>o) defined as the expectation value of the Hamiltonian 
density on \Q) in presence of the scalar condensate (pQ. We shall give an explicit formula 
for ^^(^o) which is similar to the usual perturbative expansion of the effective potential by 
means of the Feynman vacuum diagrams. Moreover, we shall show that the variational- 
perturbation theory developed in this paper offers a solution to the ultraviolet divergences 
problem in the variational approaches which is analogous to the usual perturbative renor- 
malization theory. For the sake of simplicity, we perform the explicit calculations in the case 
of selfinteracting scalar fields in one and two spatial dimensions. Indeed, these theories are 
super-renormalizable, so that we only need to renormalize the mass. In the second part of 
the paper we discuss the finite temperature corrections to the generalized Gaussian effective 
potential. Moreover we critically compare our approach to the finite temperature effective 
potential and the Gaussian potential. 
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The plane of the paper is as follow. In Section 2 we discuss the Gaussian approximation 
in scalar field theories and introduce the Gaussian effective potential. In Section 3 we set 
up the variational basis starting from the trial Gaussian vacuum wavefunctional. Section 4 
is devoted to the perturbation theory with the variational basis. The generalized Gaussian 
effective potential is discussed in Section 5. The calculations of the second order corrections 
to the Gaussian effective potential are presented in Section 6 where we discuss in detail the 
mass renormalization. In Section 7 we introduce the finite temperature generalized Gaussian 
effective potential and evaluate the lowest order thermal corrections. The second order 
thermal corrections are explicitly evaluated in Section 8. Our conclusions are draw in Section 
9. Several technical details are relegated in two Appendices. In Appendix A we perform the 
high-temperature expansions which are relevant for the lowest order thermal corrections. In 
Appendix B we collect some well-known result on the thermodynamic perturbation theory in 
the Matsubara's scheme. Moreover we present some useful results on the thermal propagator. 



In this Section we discuss the Gaussian approximation in scalar quantum field theories. 



fields in d = u + 1 space-time dimensions. 

The Gaussian approximation in quantum field theories has been widely developed since 
long time [§-§1. The Gaussian approximation is a variational method in which one considers 
trial Gaussian wavefunctionals as the ground state of the theory. 

Let us consider a real scalar field (p{x) whose Hamiltonian is 



In the Schrodinger representation the physical states are wavefunctional of (p; the conjugate 
momentum n(a;) acts as functional derivative 



II. THE GAUSSIAN APPROXIMATION 



In particular we shall focus on the Gaussian effective potential |I| for selfinteracting scalar 




(2.1) 
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The inner product is defined by 

< ^l|^2 >= / M ^2[0] . (2.3) 



The stationary Schrodinger equation reads 



<ij[(j)] = E^[(g. (2.4) 



The analogy with ordinary quantum mechanics is evident. In particular, we would like to 
apply the variational principle which has been successfully developed in quantum mechanics. 

In quantum field theories it is the ground state that determines the physical properties 
of the quantum system. The Gaussian approximation amounts to approximate the vacuum 
functional with a set of trial Gaussian functionals centered at (po: 



(2.5) 



vI/oM = exp [-- I d-'x (Fy [(t>{x) - cj),] G(f , y) [ct){y) - (j),] 

where the normalization constant is such that: 

(*o|^o) = l. (2.6) 

A well-known method to investigate the structure of a quantum field theory is to use the 
effective potential ye//(0o) 01 • In scalar field theories it turns out that the effective potential 
is the expectation value of the Hamiltonian density in a certain state wherein the expectation 
value of the scalar field is 00 0- These considerations suggested to introduce the so-called 
Gaussian effective potential Vg'((/>o)- 

The Gaussian effective potential is defined by minimizing the Hamiltonian density on 
the set of wavefunctionals Eq.(|2.5|): 



Vgep{4>o) = ^ min(^'o|i^|*o) (2.7) 

V \^o) 

where V is the spatial volume. 

VcEpiM, being a variational quantity, not only goes beyond the perturbation theory, 
but often gives a more realistic picture of the qualitative physics than the effective potential. 



Moreover the Gaussian effective potential is easily computable. To see this, we note that 
due to the translation invariance of the vacuum we have: 



G{x, y) 



Let us consider the following functional 



Afexp 



-- J (F X (F y'r]{x) G{x — y) r]{y) + 2 _/ d!" x rj{x) J {x) 



where 



(2.8) 



(2.9) 



?7(f) = 0(f) - 



and N is fixed by Eq. ( p.6|) . We can easily evaluate the following functional 



I[J] = i^il^i) . 



(2.10) 



(2.11) 



Indeed, Eq. (|2.11| ) involves a straightforward Gaussian functional integration. We get 
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I [J] = exp 



-J (Fx(FyJ{x)G-\x-y)J{y) 



2.12) 



where 



G-\x,y) 



(Fk 1 

W2g{k) 



(2.13) 



Now, we have: 



{^o\r]{xi)...7]{xn)\^o) 



(5"/[J] 



6J{xi)...6J{x'n) 



(2.14) 



J=0 



Equations p.l4| ) and ( |2.12| ) allow to evaluate the expectation values of monomial in r]{x) 
on the Gaussian vacuum functionals. It is now a straightforward exercise to calculate 



EM^,g{k)] 



(*o|*o) 



(2.15) 



We get 
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^o[</>o, m] =v{\j ^jik) + ^cpi^4 + ^ 



+ 



3 ^ 

4 4! 



d'^A; 1 



(2.16) 



(27r)-^(^)_ 

The Gaussian effective potential is obtained by minimizing i?o[0Oi fi'(^)] with respect to g{k). 
By imposing the extremum condition 

5E^[ct>,,g{k)] 



Sg{k) 







(2.17) 



we obtain 



g{k) = ^WVi^) , 



(2.18) 



where satisfies the gap-equation: 



(2.19) 



Inserting Eqs. (|2.18|) and ( p.l9| ) into Eq. ( |2.16D we get the Gaussian effective potentiah 



Vr 



GEPWO) 



^00 + ^00+27 (2vr 



'^^^ - 32 



d^k 1 



Putting 11 



we rewrite Eq. (|2.2CI| ) in the more compact form: 



(2.20) 



(2.21) 



(2.22) 



while the gap equation becomes 



2 2 ^^^2 ^ T / \ 



(2.23) 



For later convenience, it is useful to work in units of /xq = /^(0o = 0). Thus we introduce the 
following dimensionless parameters: 
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(2.24) 



A 



A 



(2.25) 





(2.26) 



Moreover we redefine the zero of the energy scale by subtracting in Vgep{4'o) the (divergent) 
quantity Vgep{4>o = 0): 



In the previous Section we have introduced the Gaussian effective potential. The most 
serious problem of the Gaussian effective potential resides in the lack of control on the 
variational approximation. For these reasons it is desirable to deal with a generalization of 
the Gaussian effective potential which allows to compute in a systematic way the corrections 
to the Gaussian approximation The problem we are interested in is not an academic 
one. Indeed, it is well known that the Gaussian approximation does not keep into account 
all the two loop contributions. As a consequence in non abelian gauge theories the Gaussian 
approximation breaks the gauge invariance 0. 

In order to evaluate the corrections to the variational Gaussian approximation we need 
to set up a variational-perturbation theory. To this end we, now, construct a variational 
basis starting from the vacuum wavefunctional ^^oM- To do this we consider "^oiv] the 
ground state wavefunctional of a suitable Hamiltonian. 

Let us consider the following operators: 




(2.27) 



III. THE VARIATIONAL BASIS 




(3.1) 
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Let |0) be the vacuum of Hq in the abstract ket formahsm. From Eq. ( |3.4| ) it follows 

(^|0) = , (3.12) 

that is ^&oW is the vacuum of Hq in the Schrodinger representation. Starting from |0) we 
can set up the many particle states by acting on the vacuum with the creation operators 
a^{p)- In the Schrodinger representation we have for instance 

r fFr p^P'^ / /■ 1 f) \ 

^iW = m = {V\a\p)\0) = / ; / d''y7^G{x,y)v{y) - 7-7^ Wo M 

J (27r)2 J2g{p} V 2 Ori[x} J 



1 1^ ^^^-y^ ^°[^] ■ (3.13) 

J (27r)2 



Obviously we also have 



Ho\i^ = [Eo + g{p}m , (3.14) 



and 



{pi\p2) = Sipi-p2){0\0) . (3.15) 

In this way we construct the orthonormal set (Fock basis) of wavefunctionals {\l'„[r7]}, where 
\l/„[?7] is obtained by applying n times the creation operator on ^I/oM- 

It should be emphasized that the Fock basis is univocally determined by the vacuum 
functional ^I/o[??]- As we will discuss in the next Section, the vacuum functional will be fixed 
with a variational procedure. For this reason the Fock basis {"^niv)} be referred to as 
a variational basis. 

IV. PERTURBATION THEORY WITH THE VARIATIONAL BASIS 

In this Section we use the variational basis to set up a perturbation theory for the ground 
state energy. To this end we split our Hamiltonian Eq. ( p.l| ) into two pieces 



H = Ho + Hi, (4.1) 



where Hq will be the free Hamiltonian and Hj the perturbation. We define Hq and Hj as 
follows: 



where 



{Ho)nm = {n\Ho\m) = SnmHnn , (4.2) 
{Hj)nm = {n\Hj\m) = (1 - Snm)Hnm , (4.3) 

Hnm={n\H\m)= j [dT]] H ^mW . (4.4) 



Equations ( [4.2|) and (|4.3| ) show that the perturbation Hi is given by the off-diagonal elements 
of the full Hamiltonian H with respect to the variational basis. If the wavefunctional 
is close to the true ground state of then we expect that the {Hj)mn are small with respect 
to {Ho)nm, i-e. Hj is a genuine perturbation. 

We recall that in Section 2 we fixed the wavefunctional W by minimizing E[(j)Q, g{k)] = 
Hqo on the class of trial Gaussian functionals, Eq. ( |2.5| ). In this way we get an optimized 
perturbation expansion. Moreover, we stress that in our scheme it is unnecessary to start 
with a small parameter in H. Thus our method goes beyond the usual perturbation theory. 

We now address ourselves in the determination of Ho and Hj. We evaluate, firstly, the 
diagonal elements of H with respect to the Fock variational basis [?]]}. To this end we 
rewrite the Hamiltonian ( |2.1|) in terms of the annihilation and creation operators ( |3.7|) and 
(p.8|). A rather length but otherwise straightforward calculation shows that: 



H = i/W + //(I) + + //(3) + i/{4) (4.5) 

where 

H^'^ = E[<f>o,g{k)], (4.6) 



2. ^.3 Xrd'^k 
m (po + -00 + -jvo 



6"" 4"^y {27rr g(k)\ 
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d^'x : ri{x) : , (4.7) 



gip) + 



+ - 



(Tk 1 



-9{P) + 



g{p) A J i2nyg{f)g{k) 

. 2 \ \j.2 I A r d"p' 1 



9{P) 



{a\p}a{-p} + a{p)a{-p)) , (4.8) 



where the colons mean normal ordering with respect to the vacuum \E'o[^]- 
In the ket formalism the n-particle wavefunctionals are given by 



\n 



(Z/J)2 

with J2i i^i = n. A straightforward calculation gives 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



(n|ifW|)n = {n\H^^^\)n = 



(4.13) 



9{Pi) + 



9{Pi) 



X r d^k 1 ^ 1 



i2nr g^k)^ 9m 



(4.14) 



As concern if'-^-*, a rather length but elementary calculation shows that in the thermody- 
namic limit V oo, we also have 



(n 



\H^*\n) = . 



(4.15) 



Now if we select the variational basis by minimizing Hqq = EQ[(j)Q, g{k)], i.e we impose the 
extremum condition Eq. ( p.l7|) we find that H^'^^ reduces to 



if*-^-* = J p g{p) a) {p) a{p) . 



(4.16) 



Moreover 
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(n 



H\n) = E[(po, 9{k)] 9{Pi) , 



(4.17) 



where g{k) = yk"^ + /i^(0o), and /i(0o) satisfies the gap equation Eq. ( |2.19| ). 

Equation ( ^.161 ) tells us that H^'^'> is the normal ordered Hamiltonian of a free scalar 
field with mass /i(0o)- Moreover it is now clear that the off-diagonal elements of the full 
Hamiltonian are due to H^^\ H^^^ and H^^\ As a consequence we can write (using again 
the gap equation): 



(4.18) 



Hi 



(Fx 



(4.19) 



We stress once again that our perturbation is given by the off-diagonal elements of the 
full Hamiltonian H with respect to the variational basis {\E'„[r;]}. This means that the 
perturbation expansion that we will discuss in the next Section, is not a weak coupling 
expansion. In other words, our variational procedure, which selects the Fock basis {\E'„[?7]}, 
minimizes the off-diagonal elements Hnm] so that even though the quartic selfcoupling is 
strong the perturbative expansion gives sensible results. Finally, it is worth mentioning that 
the simple results Eq. (|4.19|) for the perturbation Hamiltonian relies on Eq. ([4.15|) which is 
valid only for quantum systems with an infinite number of degree of freedom. 



V. THE GENERALIZED GAUSSIAN EFFECTIVE POTENTIAL 



In the previous Section we was able to split the Hamiltonian H into two pieces: the free 
Hamiltonian Hq and the interaction Hj. If we neglect Hj we see that the ground state of 
Hq is the wavefunctional "^olv] and the Gaussian effective potential Eq. ( p.22| ) is the ground 
state energy density. In other words, Vgep{4>o) is the lowest order term of the vacuum 
energy density in the perturbation expansion generated by Hj. Thus the corrections to the 
Gaussian effective potential can be readily obtained by means of the standard perturbation 
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expansion for the ground state energy. For the ground state energy we may use the well- 



known Brueckner-Goldstone formula ||1CI|| : 



EGs{<Po) = Eo{(Po,9{k)) + j: 



n=0 



mi 



Eq — Hq 



Hi |0) 



(5.1) 



For instance, up to the second order in Hi, and using Eqs. ( ^T^ ) and ([4. 3D, we have: 
Egs{(Po) =Eo + E(^oo - Hnn)-'\{n\H\0)\' . 

n>0 



(5.2) 



Higher order terms can be analyzed by means of the so-called Goldstone diagrams ||lO| , [Tl |. 



However, in order to show that Eqs in Eq- ( ^-11) gives correctly the correction to the 
Gaussian effective potential, we must ascertain that it exists a state \Q) such that: 



E, 



GSWO) 



< n 


H 


n > 


< n 


n > 



(5.3) 



with the constraint 



. 



(5.4) 



To do this, we use the Gell-Mann-Low theorem on the ground state []T2[. Let us consider 
the Hamiltonian 

H, = Ho + Hie-'\'\ e^0+. (5.5) 
Next, we introduce the temporal evolution operator 

°° ( ft H 

Ueit,to) = E —r / ^^1- / rf^ne-^(l*^l+--*-l*"l) T[Hiit^)...HiiQ] (5.6) 
nl Jto Jto 

where Hi{t) is the perturbation Hamiltonian in the interaction representation. The Gell- 
Mann-Low theorem says that if the following quantity exists to all order in perturbation 
theory: 

t/,(o,-oo)|o) _ \n) 



lim 

6^0+ (0|[/e(0,-cx))|0) {0\Q) 



(5.7) 



then 
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n > 


< 


Q > 





n > 


< 


Q > 



H ' X = ^ X • (5-8) 



Note that the denominator in Eqs. ( |5.7| ) and ( ^.8|) is crucial, for the numerator and the 
denominator do not separately exist as e 0^. 
From Eq. ( |5.8| ) it follows 

^ - -imT ■ ^ ^ ^ 

Now we show that Eqs = E. Indeed, from Eqs. ( |5.8|) and ( ^.5|) we get 

where we have taken into account that i^o|0) = -E-olO). 
Now a standard manipulation shows that 



(0|//,|fi) = m) E / ^^1- / din e[-^(*^+-+*")] X 

n=0 J — oo J~oo 

{0\T{Hj{0)Hj{t,)...Hj{Qm,onn (5.11) 

where the subscript means that we need to take into account only the connected terms. 
In order to carry out the time integrations, we consider the nth-order contribution in Eq. 
( p.ll| ). Observing that 



Hi{t) = e'^°'Hie-'"^\ (5.12) 
and that the n\ possible time orderings give identical contributions, we get 

J — oo J — oo 

By changing variables to relative times: 

Xi = ti, X2 = t2 — ti, ... ,Xn = tn — tn-1 (5-14) 

one finally obtains 
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(E - i?o)(") = (0 I Hj- ^——Hj...Hj- }^—-Hj I 0)_ . (5.15) 

Because of the limitation to connected contributions, the hmit e ^ 0+ is harmless. Hence 

we get 

„=o \Eo-Ho ) 
which shows that indeed E = Eqs- 

We can finally write down the generalization of the Gaussian effective potential we are 
looking for 



1 < VL 


H 


n > 


V <vl 


Q > 



with the constraint 



{n\n) ^^-^^^ 

Note that Eq. ( |5.18| ) assures that the expectation value of the scalar field (j){x) on the state 
\n) is 00- 

Several remarks are in order. Equation (|5.16| ) shows that E reduces to Eq in the zeTo-th 



order due to the normal ordering of the interaction Hamiltonian. Thus, in that approxima- 
tion Vg(0o) coincides with the Gaussian effective potential. 

Higher order contributions to the generalized Gaussian effective potential 1^(00) can be 



evaluated by the Brueckner-Goldstone formula Eq. ( |5.16|) . In this case one deals with an 
expansion in terms of the Goldstone diagrams [|lT|. However, one can do better if one uses 
Eqs. ( p30D and ([5Jl| ): 



Vg(0o) = E / dtn- rft„(0|T(//,(0)i//(ti)...i^/(tn))|0)_„ (5.19) 



where we have performed the harmless limit e 0"^. Indeed a given term in Eq. ( p.l9p can 



be easily evaluated by means of the standard Feynman diagrams. It is evident that the free 
Feynman propagator coincides with the propagator of a scalar field with mass /i(0o)- 
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It should be stressed that it is convenient to analyze the expansion ( |5.19|) by means of 
the Feynman diagrams. Indeed, it is well known that a Feynman diagram with k vertices 
contains k\ Goldstone diagrams, corresponding to the number of permutations of the times 
ti, tfc. 

The expansion Eq. ( p. 191 ) gives rise to a set of vacuum Feynman diagrams. Let us 
consider the term of order n in Eq. ( ^.19| ). In the time-ordered product there are n + 1 
interaction Hamiltonian. This means that, unlike the usual vacuum diagrams the factorial 
factor n\ is cancelled by the number of permutations of the independent integration variables. 
As a consequence, Eq. ( |5.19|) allows a diagrammatic expansion of the higher order corrections 
which is amenable to a diagrammatic resummation. 

Finally, we point out that in terms of Feynman diagrams, the constraint Eq. (^.18 



sets to zero the tadpole-like diagrams that are due to the linear term in the interaction 
Hamiltonian. 



VI. CORRECTION TO THE GAUSSIAN APPROXIMATION 



In the previous Section we have introduced the generalized Gaussian effective potential 
which allows to compute in a systematic way the corrections to the Gaussian approximation. 
Presently we focus on the second order corrections. Let us consider the first non trivial term 
in the perturbative expansion Eq. ( ^.19|) . We have 



with 



j"jt{0\THjiO)Hjit)\0)^ 



Taking into account Eqs. ( p.l8| ) and ( |4.19| ) we get: 



V 



dt (fx(fy 



(^)^ < I T(: rf{x) :: rf{y) :) | > + 
+ (A)'<0|r(:774(a;)::77'(y) :) | > 



(6.1) 



(6.2) 



(6.3) 
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where x = (0,x) and y = {t,y). Therefore we have (see Fig. 1): 



where the Feynman propagator is: 



(6.4) 



GF{x,y)= I (6.5) 



Inserting Eq. (|6.5| ) into Eq. ( |6.4|) and performing the time and spatial integrations we recast 
Eq. dJ) into: 



Av^c(0o) = - ^ / n 



3 



\{2ny2g{h) g[k^) + g^k^) + g{k^) 

4! J M {27vr2g{k,) g{h) + g{k,) + g{h) + 9{h) ' 
Note that the lowest order contribution in the loop expansion is the two-loop diagram 
(diagram (a) in Fig. 1) which was lost in the Gaussian approximation. In non-abelian gauge 
theories this diagram is crucial in order to maintain the gauge invariance in the variational 
Gaussian approximation 

Now we discuss the second order corrections in the case of scalar fields in = 1,2 spatial 
dimensions 1 1^ . 



A. SCALAR FIELDS IN 1+1 DIMENSION 

To start with, we consider the Gaussian effective potential in one spatial dimension u = 1: 
Vgep{<Po) = ^<^o + ^ + h{^^) - ^/o^(/i) (6.7) 



where 



^ , , 1 /■+°° dk 1 , , 

W = ^ / 7r-Y= , (6.8) 

2J-00 27r /^(^) 



1 r+oo fjU I 

IM = ^ ^V^. (6.9) 
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Introducing an ultraviolet cutoff A, it is not difficult to see that the integrals ( |6.8|) and ( |6.9|) 
display quadratic and logarithmic divergences. Subtracting the energy density of the (po = 
vacuum Vgep{4>o = 0), one is left with a logarithmic divergence which can be cured by 
renormalizing the mass. In the Gaussian approximation the renormalized mass is defined as 



2 5VGi?p(0o 

mo = 



= + ^loifio) = f^l (6.10) 

00=0 



where we recall that /io = /i(0o = 0). However we shall see that once we consider the 
corrections to the Gaussian approximation the prescription Eq. ( |6.10| ) needs modification. 



A serious problem of the variational approximation in quantum field theories is due to 



the presence of ultraviolet divergences |]14[. The variational-perturbation theory discussed in 
this paper offers a natural solution to the ultraviolet divergences problem which parallels the 
perturbative renormalization theory. As a matter of fact, we showed that the generalized 
Gaussian effective potential Vg{(Po) is the energy density of the vacuum \Q) with scalar 
condensate 0o = Observing that only the energy differences are of importance, our 

renormalization prescription will be to reabsorb the ultraviolet divergences in the field theory 
without scalar condensate. Moreover for 0o = the Hamiltonian Eqs. ( 4. IS ) and (f4.19 ) 



reduces to the one of a scalar field with mass and normal ordered quartic selfinteraction. 
In the case of one and two spatial dimensions that theory is super-renormalizable and we 
only need to renormalize the mass. We define the renormalized mass by means of the 
zero-momentum 2-point proper vertex T^'^\p; (pQ = 0): 

m^ = -r(2)(O;0o = O) . (6.11) 

In the Gaussian approximation the (po = Hamiltonian coincides with the free Hamiltonian 
of a scalar field with mass jjQ. So that Eq. (6.11) gives 

ml = fil (6.12) 



which agrees with Eq. ( |6.10| ). In one spatial dimension Eq. ( |6.12| ) eliminates completely 



the ultraviolet divergences, for the higher order corrections are finite. 



From Equation (|6.12| ) and the gap equation we get 

Ao 



m 



m. 



-loimn) . 



(6.13) 



Inserting Eq. (§J^) into Eq. (^) and using Eqs. ( ^12^^:271) and (|2]2|) one obtains |1[|: 

1 



gepWo) 



-2A$^ + 



X 



24A 



3A X — 1 

TT Z 



(6.14) 



with the gap equation 



X - 1 + - Inx = 12A$o 

TT 



(6.15) 



These results have been obtained for the first time by S.J. Chang [0. In Figure 2 display 
^GEpi^o) as a function of $o for various values of the dimensionless coupling A. Several 
features are worth mentioning. Firstly, $o is always a local minimum of Vgep{^q)- For 
A < Ac, with Ac — 2.5527, the $o = vacuum is the true ground state. On the other 
hand, for A > Ac the ground state is for $o 7^ 0. As a consequence, a first order phase 
transition occurs at Ac. However, Chang ||T6| pointed out that the Simon-Griffiths theorem 
p!7| rules out the possibility of a first-order phase transition in the one-dimensional X(f)^ field 
theory. Moreover, Chang showed that there is no contradiction between the existence 
of a second order transition and the Simon-Griffiths theorem. 

Remarkably, it turns out that the two loop correction, diagram (a) in Fig. 1, gives rise 
to a second order phase transition. Indeed that correction is finite: 

A^ 



with 



dx dy 



x 



1 

x2 + l)(y2 + l)[(a; + y)2 + l] 
1 



(6.16) 



[Vx^+1 + ViFTi + ^(x + y)2 + 1] 

So that in this approximation we get 



0.7136 



(6.17) 
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-2 - ^GEPy'^o) - O'—'^O ■ 

Lin JL 



(6.18) 



In figure 3 we display Eq. (|6.18| ). We see that now there is a second order phase transition 



at 



V2a 



0.8371. This is confirmed by considering the mass-gap of the $o = vacuum: 



m: 



phys 



^0 + S(0) 



(6.19) 



where S(p) is the proper self-mass of the 0o = theory. In the second order approximation 
S(p) is given by the so-called setting sun diagram |T8[. It is easy to show that 

A2 



j:{0) = --Jd^x[GE{x)f 



(6.20) 



(Pk e^^^ 
(27r)2 P + 



where Ge{x) is the Euclidean Feynman propagator 

Ge{x) -- 

A straightforward calculation gives 

E(0) = -2aAVo • 
Thus, from Eqs. ( |6.19|) and ( |6.22| ) we get 



(6.21) 



(6.22) 



phys 



, AL 



(6.23) 



One can easily check that in this case 



m: 



GWO) 



phys 



(6.24) 



<l>o=0 



A remarkable consequence of Eq. (|6.24|) is that the mass renormalization of the Gaussian 
effective potential extends to Vg{4'o) in the two-loop approximation. Equation ( |6.23| ) tells 
us that the (pQ = vacuum is stable for A < Ac. Moreover near the critical coupling we have 



^^phys 



(Ae-A)2 



(6.25) 



so that the correlation length ^ = — - — diverges as 

'^phys 
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e ~ (Ae - A)-^ v=\. (6.26) 



Our results are in agreement with previous studies [p!6| , p!9| , p0| . However our generalized 
Gaussian effective potential relies on a firm field-theoretical basis which allow us to take 
care of the higher order corrections. Moreover in our scheme there are not ambiguities in 
the renormalization of the ultraviolet divergences. 

Let us consider now the contribution due to the diagram (b) in Fig. 1. We have 

AyG(0o) = -a — ct)l^ll - h— (6.27) 

X X 



where 



3 1 



Vx^+l,/^i]^Tly/^+lJ{x + y + zy + l 



0.0509 . (6.28) 



We would like to stress that now Eq. ( 6.24 ) is no longer valid. In the present case this does 
not matter due to the fact that the second order corrections are ultraviolet finite. However, 
in the case of two spatial dimensions these corrections are divergent. Thus, adopting the 
renormalization prescription Eq. (|6.24 ) instead of Eq. ( p. 11 ) it may lead to an incongruous 
result. 

In Figure 4 we contrast the generalized Gaussian effective potential in the two-loop 
approximation (full lines) and in the full second order approximation (dashed lines). 

Few comments are in order. The order of the transition is not modified by the second 
order three- loop correction. Moreover the critical coupling Ac — 1.1486 is quite close to our 
previous value. As a matter of fact, in the critical region A ~ 1 the effects of the three- 
loop correction do not substantially change the shape of the potential. Therefore we can 
safely conclude that the most important contributions in the critical region are given by the 
two loop term. This suggests that higher order corrections do not modify the order of the 
transition. 
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B. SCALAR FIELDS IN 2+1 DIMENSIONS 

In the case of two spatial dimensions the Gaussian effective potential is 

W(0o) = + + - ^/o'(^) , (6.29) 



where, now, 



IM = t; / 7^^^ (6-30) 



Subtracting the energy density of the 4>o = vacuum we get 

W(0o) - W(0) = + ^c^t + A(^) - /i(^o) - ^[/o'(^) - ^'(/^o)] , (6.32) 

which is still affected by ultraviolet divergences. 
Introducing the renormalized mass 



2 d'^VGEp{(i)o) 

m, — 



A 



= + -/o(^o) = Ato (6-33) 

00=0 



and using the gap equation 



f^\M = m' + ^cl^l + ^l'M (6.34) 



we get the finite result M 



"^0 , ^^4 (v^-1^^ ' 



2 " 2A-K 

Similarly we can rewrite the gap equation Eq. ( |6.34| ) as 

3 



l + ;^A + 2v^ 
27r 



whose explicit positive solution is: 



3r^ 



(6.35) 



a;= l + 12A<l'g--A(v^-l) (6.36) 

TT 
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In Figure 5 we show Vj^p($o) versus $o for three different values of A. Again we find 
a first order phase transition at the critical coupling A ~ 3.0784 Note that in two 
spatial dimensions there are not rigorous results which could exclude a first order transition. 
Nevertheless it is important to investigate the effects of the second order corrections to the 



Gaussian effective potential. From Equation ( |6.4| ) we have 



AVg{%) = -I dx J d'z + (6.38) 

where z = {t,x). Now, observing that the Feynman propagator is an even function and 
performing the Wick rotation, we obtain: 

^Vg{%) = -\j d'zE I^GUze) + ^GUze)] (6.39) 



where 2 1 



with z = \ze\ and Ki is the modified Bessel function of order 

i^.(x) = yje-. (6.41) 
Thus we get for the Euclidean Feynman propagator: 

-flZ 

Ge{ze) = —- . (6.42) 
At:z 



Introducing 



/3(/i) = / d^z, GUze) = / (6-43) 



hif^) = I d'zEG%{zE) = I , (6.44) 



we rewrite Eq. (|6.39|) as 



A^g(0o) = -^^3(^) - ^/4(^) . (6.45) 
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From Eqs. (|6.43| ) and ( |6.44| ) we see that J3 and I4 are divergent. We regularize the integrals 



by means of an the ultraviolet cutoff e 
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1 r°° , e 
dz — 



167r2 



167r2 



[ln(/ie) +ln3 + 7] + 0(e) 



(6.46) 



(4^) 



00 

dz — 



{Att 



^ ^i + 4^[ln(^e)+ln4 + 7-l][> + C>(e) 



(6.47) 



where 7 is the Euler-Mascheroni constant. 
Putting it all together we obtain 

Vg{<Po) - Vg{0) = + + IM - Ji/xo) - ^[PM - + 



+ ^[ln(/.e) + ln3 + 7] 



1927r2 
' 7687r3 



7^(^-^o)[ln4 + 7-l] + 



[// ln(/ie)] - //oln(/ioe)] 



(6.48) 



Now we show that the logarithmic divergences are cured by renormalizing the mass. To this 
end, we observe that 111 



(6.49) 



Inserting Eqs. ( |6.49| ) and ( |6.50| ) into Eq. ( |6.48| ), and using the gap equation Eq. 
rewrite Eq. ( |6.48[ ) as: 

1 



^g(0o) - ^g(O) = + 



/^o 

2477 



A 



1287r2^°» \/ ,.2 



1927r^ 



[ln(/ie) + In 3 + 7] 



A^ 



- //o)(ln4 + 7 - 1) + /i ln(/ie) - /io ln(/ioe)] 



7687r3 

As we have already discussed, the renormalized mass is 



(6.50) 



we 



(6.51) 



m^ = -r(0,(/)o = 0) = /i^ + S(0) 



(6.52) 
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In the lowest order Gaussian approximation we have S(0) =0 and Eq. ( |6.52| ) reduces to 
Eq. (|6.33| ). In the second order approximation we must introduce a mass counterterm so 
that 

S(0) = W - ^ 1 d'xE GUxe) (6.53) 

where the second order term is due to the setting-sun diagram. Explicitly, by using the 
previous regularization, we find 

S(0) = W + ^[ln(^oe) + ln3 + 7] + 0(e) ■ (6.54) 

We fix the mass counterterm by imposing that 

ml = fil. (6.55) 

This results in: 

= + ln3 + 7] . (6.56) 

As a consequence, we must introduce the following counterterm Hamiltonian in the (po = 
theory: 

SH=^ 6m^ J d^x(f)\x) . (6.57) 

After writing (p{x) = (po + rj{x) and using the constraint Eq.( ^.18| ), it turns out that 5H 
adds to the second order generalized Gaussian effective potential the further contributions 
depicted in Fig. 6: 

5V^^^ = ]-(t)l 6m^ (6.58) 



5V§^ = lsm^f . (6.59) 



Now, observing that 



2 J i2n)^2g{k) 



(6.60) 



2g{k) J-^ 27r kl + g'^(k) ' 
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we have (e — > 0): 

1 



^2 ^ ^ Stt 



+ 0(e) . (6.61) 



Finally, using Eq. ( |6.56|) we get: 



SV^''^ = -^[ln(/ioe) + ln3 + 7] (6.62) 



Aire 



(6.63) 



Now, it is easy to see that ^V^."^ eliminates the ultraviolet divergences of the two-loop second 
order correction, whereas 6Vq^ makes finite the three-loop second order correction. Thus 
we are left with the finite result: 



2 



+^lZLln +^l-^[l-ln(i^)-ln-] 
A^ 4 

r/^o[l-ln-]. (6.64) 



In terms of the scaled variables Eqs. (|2.23| ) and (|2.24|) we rewrite Eq. (|6.64| ) as: 



Voi^o) = V^tU^o) + In + ^(v^ - 1)(1 - In -) - ^A^V^ In • (6.65) 



It is worthwhile to study separately the effects of the two second order corrections. If we 
take into account the two-loop correction we get 



l4'^)($o) = VS^U'^o) + ^$^lny^ . (6.66) 



3A 



2 



7r2 



In Figure 7 we display Eq. ( |6.66| ) for A = 1, 3, and 5. As it is evident there is no spontaneous 
symmetry breaking |jl3|]. Comparing Fig. 7 with Fig. 5, we see that the two-loop correction 
adds to the Gaussian effective potential a positive contribution which is important in the 
region $o ~ 1 and overcomes the negative minimum displayed by Vq^p{^o) in that region. 
On the other hand considering the three-loop second order correction, we have 
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vi'\<l>o) = V^^W^o) + ^(v^- 1)(1 - In^) - ^A^v^lnv^ . (6.67) 

From Fig. 8, where display Eq. (|6.67|) for three different values of A, we deduce that the 
most important effects of the three-loop second order correction is near the origin where one 
gets: 

Indeed, Fig. 8 shows that Vq\^o) undergoes a continuous phase transition at Ac = 3.0959. 
This feature persists even for the full second order generalized Gaussian effective potential 
Eq. ( |6.65| ) (see Fig. 9). Our result is in qualitative agreement with Refs. |^ and pl| . 
However, from Fig. 9 we see that the condensation energy is very small. Moreover the 
critical coupling Ac = 3.0959 differs from that of the Gaussian effective potential by less 
than one percent. We feel that the only sound conclusion we can draw is the exclusion of a 



first-order phase transition. Note that unlike of what stated in Ref. [|J], the absence of a 



broken phase is not in contradiction with the analysis by S.F. Magruder and S. Chang 



and S.F. Magruder |26|. 

We would like to conclude this rather technical Section by stressing the most important 
achievements of our analysis. Our analysis of the ultraviolet divergences in two spatial di- 
mensions showed that our renormalization procedure works up to the second order. However, 
it is clear that our renormalization can be extended to the higher orders by the usual renor- 
malization procedure. Thus we feel that our results put the generalized Gaussian effective 
potential on the same level as the effective potential. 



VII. THERMAL CORRECTIONS TO THE GENERALIZED GAUSSIAN 

EFFECTIVE POTENTIAL 

The aim of this Section is to study the thermal corrections to the generalized Gaussian 
effective potential. For reader convenience, let us firstly recall the essential points of the 
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finite temperature effective potential p7|-pg[] and tlie finite temperature Gaussian effective 
potential MM]. 



Following the classical paper by L. Dolan and R. Jackiw the finite temperature 
effective potential in the one-loop approximation is given by: 



1 r d^k 



HE' + ^l) 



(7.1) 



where Un = 2nf3n, /3 = -k, are the Matsubara's frequencies and E'^ = k"^ + M 



)) = + f 00- Performing the sum over n one finds p9| : 



r- 



"(0o) 



d^k E 1 



d^k 



ln[l 



-PET. 



(7.2) 



^ J {27iy 2 ' pj {27iy 

In the right hand of Eq. ( [7.2|) the first term is the zero-temperature one-loop effective 
potential, while the second term gives the one-loop thermal corrections. 

As concern the Gaussian effective potential at finite temperature, we shall follow G.A. 
Hajj and P.M. Stevenson Let us consider a system in thermal equilibrium; this means 
that our system has minimized its free energy: 



where Z is the partition function: 



(7.3) 



(7.4) 



In order to evaluate the thermal corrections to the Gaussian effective potential, the authors 
of Ref. 1^ split the Hamiltonian as 



H = Ho + Hj 



(7.5) 



where Hq is the Hamiltonian of a scalar field with variational mass M^, while Hj comprises 



the remainder. The variational mass is fixed by minimizing the free energy Eq. ( |7.3| ). To 
do this one uses the thermodynamic perturbation theory to evaluate the free energy in the 
lowest order in the perturbation Hamiltonian. Writing 
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we get 

Z~Tre-^^°[l-/3(i//)^] (7.7) 

where (O)^ means the thermal average with respect to Hq: 

, Tr(e-/^^°0) 

- Tr(e-/^^o) • ^^-8) 

From Eqs. ( [7.3|) and ( |7.7| ) we obtain: 

F = In Tr (e"^-^") + < Hj . (7.9) 

Now observing that the thermal average involves a summation over the eigenstates of Hq, 
it is not too difficult to find 



^gV(0o) = ^ = /i + If - ^(/o + I'o? + + (7.10) 

where 

^-^/(C^^^(^--"^^^o^ ^^•^^) 

f d'^k 1 1 
° y (27r)-^(^)e/3sW-l ^ • ^ 

with g'(fc) = \Jk'^ + . It turns out that the mass M satisfies the thermal gap-equation: 

M' = m' + ^[Io + I^ + cl^l]. (7.13) 

A remarkable consequence of Eq. ( [7.10| ) is that the finite temperature Gaussian effective 
potential can be obtained from Vgep{4>q) with the substitution rule: 

lo^Io + I^ (7.14) 



Ji^/i + /f. (7.15) 
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The main drawback of the Hajj and Stevenson's approach is that the sphtting of the Hamil- 
tonian in Eq. ( [7 .51 ) is not natural, for the variational mass M, which is fixed by minimizing 
the free energy, depends on the approximation adopted in evaluating perturbatively the free 
energy. Moreover, the calculations of the thermal corrections beyond the Gaussian approxi- 
mation is very difficult. On the other hand, as we have already discussed in Section 4, in our 
approach the Hamiltonian is split into two pieces, the free Hamiltonian and the interaction, 
in a natural manner. 

As a consequence the thermal corrections to the generalized Gaussian effective potential 
can be evaluated easily by means of the familiar thermodynamic perturbation theory |3^. In 
the remainder of this Section we focus on the lowest order thermal corrections and compare 
with the one-loop thermal effective potential corrections and the finite temperature Gaussian 
effective potential. The higher order thermal corrections will be discussed in the next Section. 



In the lowest order in the perturbation we write 

1 



In Tr(e 



-fSHo 



+ {Hi) 



(7.16) 



where, now, Hq is given by Eq. ( 4.18 ) and Hj by Eq. ( 4.19| ). Observing that the eigenstates 
of Hq are the states \n), Eq. ( |4.11|) , with eigenvalues En given Eq. ( [4.171 ), we have: 



According to our definition Eq. ([4.3|) we have {n\Hi\n) = 0, so we end with 



(7.17) 



{HiY 







(7.18) 



The calculation of the partition function Zq = Tre is straightforward: 



-l3Eo 



n 



1 _ Q-/3g{k) 



(7.19) 



where g{k) = yk'^ + ^'^{(po), fi^{(po) satisfies the zero temperature gap equation Eq. ( p. 191 ) 



and ^0 = E[cf)o, g{k)], Eq. The insertion of Eqs. ( [TTBI ) and ( frTPI) into Eq. 

leads to 
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Vg{<Po) = ^ = Vgep{<Po) + ^J^ 1^(1 - e""'^'^) ■ (7-20) 



Note that Eq. (|7.20|) differs from the finite temperature Gaussian effective potential Eq. 



( [7.16| ). The difference resides in the different use of the gap equation. In our scheme the 
gap equation Eg. ( p.l9| ) is fixed once and for alL In particular it does not depend on the 
temperature. On the other hand, in the finite temperature Gaussian effective potential 
approach the gap equation includes the thermal effects. The gap equation fixes the basis to 
sum over in the thermal average, so that different gap equations lead to inequivalent basis. 
In fact the discrepancy between our results Eq. ( |7.20| ) and the finite temperature Gaussian 
effective potential comes from the thermal average of the interaction Hamiltonian. In our 
approach Eq. ( [7. 181) holds, whereas in Ref. |^ {Hi)^ 7^ 0. Note that the possibility of 



non-equivalent basis is a peculiar feature of quantum systems with an infinite number of 
degrees of freedom. 

It is worthwhile to compare the one-loop thermal correction to the effective potential 
with our finite temperature generalized Gaussian effective potential Eq. ( [7.20|) . Comparing 
Eq. (|7.2|) with Eq. ( |7.2CI| ) we see that the former agrees with the latter if 



M\<P,)=m' + ^(g^,,\(t^o). (7.21) 

Now /i^(0o) satisfies the gap equation Eq. ( |2.19| ) which, from a diagrammatic point of view, 
is obtained by summing the infinite set of the superdaisy graphs in the zero temperature 
propagator. In other words, if in the thermal corrections of the effective potential we replace 
the tree level mass of the shifted theory with the mass /i^(0o) obtained by summing the 
superdaisy graphs at T = in the propagator, then we obtain again a free energy density. 
Up to now this remarkable result in thermal scalar field theories holds for the one-loop 
approximation. In the next Section we will show that it extends to higher order thermal 
corrections too. 

Let us analyze the lowest order thermal corrections Eq. ( [7.20D in the case v = 1,2 
In one spatial dimension Eq. ( |7.2(]| ) reads: 
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VS{%) = V^t'p + ^,rdt In [1 - e"v/^] (7.22) 

where P = /3/io- In Figure 10 we show Vq{^o) — Vq{0) (in units of fi^) versus $o for A > Ac. 
As we can see, the symmetry broken at T = gets restored for T > T^.- Obviously the 
critical temperature depends on A. For A = 4 we find Tc — 1-27. It turns out that Tc can be 
estimated, within a few percent, by means of the high-temperature expansion of the integral 
in Eq. ( |7.22|) . From the results of the Appendix A (see Eq. ([A16| )) we find the following 



high-temperature expansion: 



yJ^($o) = Vhfpi^o) - ^ + ^ + ^ln(^ 

C(3)^2a2 , C(5) .4 . , .0.^6 



where Q{z) is the Riemann's zeta function. In Figure 10 we also show the high-temperature 
expansion Eq.( |7.23| ) (dashed lines); as we can see the high-temperature expansion is a very 



good approximation even near the critical temperature. Indeed, for A = 4 Eq. (|7.23| ) predicts 
a critical temperature which differs by less than one percent from the numerically estimated 
value. 

The case of two spatial dimensions can be dealt with in a similar way. We have : 

Vj(*o) = V^fp{^o) + -^J^dtt ln(l - e-v/^) . (7.24) 

27r/33 Jo 

In Figure 11 we display Eq. ( [7.24| ) (we subtract the temperature dependent constant Vq{0)) 
for three different values of the temperature and A = 4. Again the thermal corrections lead 
to the expected symmetry restoration at high temperatures. As in the previous case we 
performed the high-temperature expansion of the integral in Eq. (|7.24]) . We find (see 
Appendix A): 



= V^tU^o) + ^ - ^ln(/3^x) + r^ln^^ + 4^ - ^^-^ . (7.25) 

However, we would like to stress that the expansion parameter in the above mentioned 
integral is f3^x. So that in the region $o ~ 1 where x ^ 1 the high temperature expansion 
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Eq. (|7.25| ) breaks down. In Appendix A we develop an alternative expansion which is useful 
in the region 13'^ x > 1. 

To conclude this Section, it is worthwhile to perform a quantitative comparison of our 
generalized Gaussian effective potential with the finite temperature Gaussian effective po- 
tential and the one-loop thermal effective potential. For definiteness we focus on the critical 
temperature as a function of the coupling constant A in the case of two spatial dimensions. 
In this case the one loop thermal effective potential reads (assuming unit mass): 

V^'^'^'m = + AV + ^(1 + 12X^1) - Y^(l - 12A$^)t - ^ + /f (1 + 12A$^) . 

(7.26) 

As concern the finite temperature Gaussian effective potential, the critical temperature can 
be extracted from Eqs. (|7.1CI| ), (|7.11| ) and (|7.12|) with u = 2. In Figure 12 we compare the 



critical temperature as a function of A (in units of Ac) for the three different potentials. We 
see that our finite temperature generalized Gaussian effective potential leads to a critical 
temperature which increases more slowly than for the other two potentials. This is due to 
our choice of the variational basis which implies (Hj)^ = . 

VIII. FINITE TEMPERATURE DIAGRAMMATIC EXPANSION 

In the previous Section we evaluated the lowest order thermal corrections by means of the 
thermodynamic perturbation theory. Presently we would like to calculate the higher order 
thermal corrections. To do this the usual thermodynamic perturbation theory is useless. 
Instead we may follow the Matsubara's methods p2| , |33| . In the Matsubara's scheme one 
deals with scalar fields which depend on the fictitious imaginary time r varying in the 
interval (0,/?). If the Hamiltonian of the system in thermal equilibrium can be written as 
H = Hq + Hj, then one can show that the corrections to the thermodynamic potential are 
given by (see Appendix B): 



Afi = In (T, exp |- jj^^ i//(r)dr| 



)' (8.1) 
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where Hi{t) is the interaction Hamiltonian in the Matsubara's interaction representation. 
Tj. is the r-ordering operator and the thermal averages are done with respect to the free 
field partition function. Moreover, it turns out that only the connected diagrams contribute 
to AQ: 

(3 ^1 ml Jo 
In our case, if we write 

V^iM = VcEPi^o) + ^ / 1^ (1 - ^"''^'0 + ^^oi^o) , (8.3) 
we readily get: 

Note that, due to Eq. ( [7.18| ), the sum in Eq. ( |8.4| ) starts from m = 2. The thermal average 
of the time-ordered products is evaluated by means of the Wick's theorem for thermal fields 
P5| . In this way we obtain the thermal corrections to the generalized Gaussian effective 
potential by means of the connected thermal vacuum diagrams. For instance, the second 
order thermal corrections are displayed in Fig. 13. The vertices can be extracted from the 
interaction Hamiltonian Eq. ([4.19|) . The solid lines in Fig. 13 are the thermal propagators 
of the free scalar fields with mass /i(0o)- 

Gp{x -y,Ti- T2) = < T^riix, Ti)r]{y, T2) >^= - / TTTT^ ; ;7r^ (8-5) 

where uin = 27rn(3. Note that the thermal propagator is periodic in the time variable with 
period 2-11(3. 

A distinguishing feature of the graphical expansion of Eq. ( ^.4| ) with respect to Eq. 
( p.l9| ) stems from the fact that the factor (m!)~^ coming from the mth order term is not 
completely cancelled by the number of different Wick contractions corresponding to a given 
graph. Consequently, a graph contributes to AVj^(0o) in proportion to a combinatoric 
coefficient depending on the order of the graph. Moreover, in evaluating the contribution due 
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to a given graph one should take care of the normal ordering prescription in the interaction 
Hamiltonian. In turns out that the normal ordering in Hi modifies the so-called anomalous 
diagrams, i.e. the diagrams which vanish at zero temperature 



For instance, in Fig. 13 



the diagrams (b), (c) and (e) are anomalous. To see this, we note that the normal ordering 
in Hi is ineffective when we consider a thermal contraction of two scalar fields belonging 
to different vertices. Therefore, the normal ordering comes into play when we contract two 
fields which belong to the same vertex. In this case we get the following thermal average: 



G^(0) = (T.:r](f,r)?7(f,r) :)^ 



.6) 



instead of Gp{Q). Taking into the account canonical commutation relations between the 
creation and annihilation operators it is straightforward to show that 

dJ'k 1 



Now we observe that 



G'^(O) = G^(0) 



(Fk 1 



lim ^^3(0) . 



(27r)- 2g{k) 
Indeed, from Eq. (|8.5|) it follows that: 

1 +°° r rl'^h 1 



f3^^^J {2'KYul + g\k) 



By using the well known identity |2T 

cotgh(a;) = 

we rewrite Eq. (|8.9|) as: 

- 



1 2x ^ 1 

TTX TT + in? 



d^k 1 



-cotgh 



Pgik) 



{^^r2gik) 

Finally, performing the limit /3 — > oo in Eq. ( ^.11| ) we recover Eq. 
Using Eq. ( ^.11|) we rewrite Eq. ( ^^) as: 



G^(0) 



d^k 1 
W2g{k) 



coth i^g{k)\ - 1 



(8.7) 



L8) 



(8.9) 



^.10) 



(8.11) 



(8.12) 
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Note that Gf^{0) is free from ultraviolet divergences in any spatial dimensions. This means 
that the ultraviolet divergences due to the tadpole Gp^O) are cured by normal ordering the 
Hamiltonian at T = 0, in accordance with the well known result that the thermal corrections 



in quantum field theories are ultraviolet finite |37 



We are, now, in the position of extending the result implied by Eq. ( [7.21| ) to the higher 



order thermal corrections. To this end we observe that the higher order thermal corrections 
to the effective potential are given by Eq. (1.9) of ref. Observing that in the imaginary 



time formalism the interaction lagrangian agrees with the interaction Hamiltonian and that 
the Gaussian functional integrations with periodic boundary conditions in Ref. corre- 
spond to the thermal Wick theorem, we obtain the desired result . There are, however, 
two further points which need to be discussed. First, our interaction Hamiltonian is normal 
ordered at T = 0. However, our previous discussion tells us that the normal ordering does 
not affect the thermal corrections, for ^^3(0) and G'/3(0) differ by a temperature indepen- 



dent term. Secondly, in Ref. there is not the linear term in the shifted scalar field. 
This means that our substitution rule Eq. ( |7.21| ) holds for the physically relevant on-shell 
thermal effective potential. 

Let us, now, explicitly evaluate the second order thermal corrections in the case of one 
spatial dimension. From Eq. ( ^.41) we have: 

1 

^V^{<Po) = I dn dT2 {T,Hj{n)Hj{T2))i^^ , (8.13) 

which gives rise to the diagrams depicted in Fig. 13. 

It is easy to see that graph (a) is temperature- independent. So it does not contribute to 
AVQ{(f)Q) due to the stability condition {Q\r]\Q) = 0. As concern the graph (b) we get: 

{b) = (^/^Vo - ^00^ J^^ dxdy dTidT2{Trr]{x,Ti)ri{y,T2))^ 

{T,:r]{x,T2My,r2)-y . (8.14) 



According to our previous discussion we obtain: 



W = (1 - ^0^) G,(0) . (8.15) 
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In a similar way we find: 



8 



For the graph (d) we have: 

id)- 



+00 



Using Eq. (|8.5|) and the result (see Appendix B, Eq. ( |B32|) ) 



-9{k)\r\ 



P^ujl + g^{k) " 2g{k) ■ 2g{k) 



+ 



\Q9{k)T _|_ ^~g{k)7 



we get 



dki dk2 



48(27r)2yo ^'7-00 g{ki)g{k2)g{k^)]\[ 



where Y^f^i ki = 0. Finally, using again Eq. ( |8.18| ) we find 



32(27r) 



Gj(0) I ' dr 



dk 



-g{k)T 



+ 



'g{k)T _|_ ^~g{k)T 
QPg{k) _ I 



(8.16) 



(8.17) 



(8.18) 



(8.19) 



(8.20) 



(/) 



384(27r)3 Jq 



dr 



+00 



dkidk2dk3 



g{.h)g{.k2)g{k^)g{ki 



X 



n 

j=i 



'g{ki)T 



og{ki)r 



+ 



+ e 



-g(ki)r 



oPg(ki) 



.21) 



with Eti h = ^- 

Some comments are in order. In Equations (8.19-21) the r-integration can be performed 
explicitly, while the remaining integrations over the momenta must be handled numerically. 
In the limit (3 — > 00 (T ^ 0) the anomalous graphs (b), (c) and (e) go exponentially to 
zero due to the factor (5/3(0). On the other hand, the graphs (d) and (f) reduce to the zero 
temperature second order generalized Gaussian effective potential. Indeed, in that limit in 
Eqs. (8.19-21) only the factors e"^*^^'-*"^ survive. Performing the elementary r-integration we 
obtain the zero temperature contributions. As a consequence the zero temperature limit of 
\/J(0o) reduces to Vg(0o)- 



37 



In the high temperature hmit we find that the graphs (e) and (f) dominate: 

Therefore, in the intermediate temperature region /3 ~ 1 we expect that the main contribute 
to Vj^(0 o) comes from the graphs (d), (e) and (f). This is indeed the shown in Fig. 

14 where we display the contributions due to the second order graphs for (3 = 1. 

In Figure 15 we display the finite temperature generalized Gaussian effective potential 
(in units of /Xq) for three different values of T and A > Ac- We see that the symmetry broken 
at T = gets restored by increasing the temperature through a continuous phase transition. 

We have also performed the analysis of the second order thermal corrections in two 
spatial dimensions. The calculation are very similar to the previous case. Moreover we find 
that the contributions to the second order thermal corrections behave similarly to the ones 
of the one- dimensional case. So we do not discuss any further this matter. 

IX. CONCLUSIONS 

In this paper we have developed a perturbation theory with a variational basis for self- 
interacting scalar quantum field theories. Our aim was to evaluate in a systematic manner 
the corrections to the variational Gaussian approximation. In particular we introduced the 
generalized Gaussian effective potential which allowed to determine the corrections to the 
Gaussian effective potential. 

Our method has been illustrated in the case of a selfinteracting scalar fields. However 
we feel that there are no problems in extending our method to scalar fields with continuous 
internal symmetry. As a matter of fact, recently our approach has been applied to scalar 



fields with 0(2) internal symmetry p9|. 



One of the most serious problem of the variational approximation in quantum field the- 
ories is due to the apparency of the ultraviolet divergences. The variational-perturbation 
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theory developed in the present paper offers a solution to the ultraviolet divergences prob- 
lem which is similar to the well known perturbative renormalization theory. Indeed, starting 
from the fact that the generalized Gaussian effective potential by definition is the vacuum 
energy density in presence of scalar condensate, we showed that the divergences are cured 
by the counterterms of the underlying field theory without scalar condensate. 

We would like to stress that, to our knowledge, there are no rigorous results on the 
problem of ultraviolet divergences in the variational approach to quantum field theories. 
For this reason we focused on scalar field theories in one and two spatial dimensions, where 
one only needs to renormalize the mass. In one spatial dimension we showed that the lowest 
order renormalization of the mass assures that the higher order corrections are finite. In the 
case of two spatial dimensions we find that the our mass renormalization procedure works 
up to the second order. However, it should be clear that our prescription can be extended 
to the higher orders without problems. 

In the second part of the paper we studied the thermal corrections to our effective 
potential. In particular in our method the Hamiltonian is split into a free piece and an 
interaction in a natural way. This allows us to directly use the well developed thermodynamic 
perturbation theory to evaluate the thermodynamic potential. A remarkable consequence of 
our analysis is that the thermal corrections to the generalized Gaussian effective potential 
agree with the ones of the effective potential provided we use Eq. (|7.21|) . 

Let us conclude by briefly discussing the more realistic case of scalar fields in three spatial 
dimensions. There is a growing evidence that quartic selfinteracting scalar field theories are 
trivial in four dimensional spacetime [^. However, recently M. Consoli and P. M. Stevenson 



proposed that the vacuum of the {Xcl)'^)^ theory is not trivial |^T|. More precisely, within the 



Gaussian variational approximation they argued that the elementary excitations behave as 
free fields while the vacuum resembles a Bose condensate. 

Recently, this triviality and spontaneous symmetry breaking scenario found some evi- 
dence in the lattice approach [^,0 • If this turns out to be the case, we expect that the 



symmetry broken at zero temperature gets restored by increasing the temperature. Thus 
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our approach to the calculation of the thermal corrections may be useful to investigate the 
nature of the thermal phase transition. In particular it is important to ascertain if the phase 
transition is first order or continuous. 



APPENDIX A: 

We are interested in the high-temperature expansion (3 —>■ oi the following integral 



h{a') 



dt In 



7r/?2 Jo 

where = . Following Ref we consider 

Bh 1 



1 - e 



(Al) 



dt 



da? 27r/32 h ./WT^{e^^+^ - 1) ' 
To perform the high-temperature expansion of Eq. (^) it is useful to deal with 

hi{e,a^)= / dt , = — , ^ ^ 

^ ' Jo VWT^{e^^-l) 



(A2) 



(A3) 



with e ^ 0^. Using the identity ^1 



— n = 1 — cothfyr-y) , 



(A4) 



we rewrite Eq. ([A3| ) as 



hi{e,a) = / 
Jo 



oo ^~ 

dt- 



VW+a^ 



E 



Vt^+a? 



t2 + a2 + 47r2n2 2 



(A5) 



Performing the change of variable y 



Va'^+iir'^n'^ 



of Eq. 



as: 



+00 



1 



^2 + 47r2n2) 

The last integral can be performed to yield ||21||: 



we rewrite the first term in the right hand 



(A6) 



dy— 

y + 1 



2 X 1 + e 
2 



1] . (AT) 



Using the definition of the Riemann's zeta function we get: 
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/«(a) = ^ + 2-i-V-^C(l + e) + /(a) = ^ + ^ + ^(^ - In 27r) + /(a) 



where 



1 °° 1 

^ n=l " 



1 + 



As concern the integral if^o?) in Eq. (^) we have [21 



1 , x~ 
dx- 



2 Jo y/x^ 



1 r(i-.)r(f) 
2" 21-^ r(i-f) 

- ^ + i In a . 



So that in the hmit e — > O"*" we obtain: 



TT 1 1 
hm /ii(e,a^) = — + - In a+ -[7 - ln47r + /(a^) 

e^0+ Zci 2 2 



Finally we perform the Taylor expansion of I{a ): 



C(3) o , 3C(5) , 



167r2 2567r4 



Putting it all together we obtain: 



TT 1 



hi{a') = — + -lna + -(7-ln47r) 



C(3) 2 , 3C(5) 
167r2 256 vr^ 



Whence: 



i^) = ^ I - + 1 Ina + -(7 - ln47r) - + + 0{a') 

27r/32\2a 2 2^' ' 16%^ 256 tt^ ^ ^ 



In order to recover /i(a2) we integrate h^{a'^) in with the boundary condition 



^(0) = ^J dt ln(l - e-*) 

7r/32 Jo 



TV 



6/?2 



We get: 



h{a^) = — <i — + vra - — + — In ( — I + 



27r/32 [ 3 



4 2 V47r 



^ C(3) + ;r3^C(5) + 0(a«^ 



327r2 



2567r4 
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Let us, now, evaluate the high-temperature expansion (a 0) of the following 



poo ^-p- p- 1 poo I — 

J{a^) = dtt ln(l - e-^^^) = - / dy ln(l - eV^+'^ ) 
Jo 2 JO 



To this end, we evaluate 



dJ 1 



1 



Using the identity Eq. (|A4|) we write 



J'(a^)= hm[K,^(a^) + ir,'^(a^)] 



where 



K. 



1 1 



E 

oo 



4 f-^ Jo y + a? + An^n"^ 



I fCO y 

dy 



VyT 



To evaluate K^^^a^) we proceed as we did for I^^^a^). We obtain 

1 



4e 



' +-4-C(2e)--|^eC(2) 



a2e (47^2)^ 

lna+|?,C(2) 



(47^2)1+. 

0(e). 



As concern K^\a?') we find 
Using (^(2) = ^ we obtain 



.l-2e 



J'(a2) = -1 In a + 



1 



a a 
4 ~ 48 



Integrating in we are led to 



,3 „A 



J{a^) = J(0) + -a^ - iaMn + — - — . 
^ ^ ^ ^ 4 4 6 96 



where 
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J(0) = / rftHn(l - e"*) = -C(3) (A26) 
Jo 

It is useful to perform also the low-temperature expansion (a oo) of J {a?). To do this we 
note that: 

1 /-oo p-n^'yW- 

Aa') = -1,1: dy . (A27) 

Changing the integration variable we get 

„2 00 1 00 —na 00 —na 

J(''')--jT.-l dte-'-^^-aj:'— - Y.— - (A28) 

^ n=l ■^'^ n=l " n=l 

This last expression can be used to approximate J(a^) for a > 1. 

APPENDIX B: 

For reader convenience we briefly discuss the thermodynamic perturbation theory in the 



Matsubara's scheme [11,35]. Let us suppose that the Hamiltonian of our thermodynamic 



system can be written as 

H = Ho + Hi . (Bl) 

We are interested in evaluating the thermodynamic potentials perturbatively in Hj. To this 
end we introduce the 5-matrix: 

e^^ = S-\T)e"°\ 0<T<p. (B2) 

Let us consider the field operators in the Matsubara's interaction representation: 

<j){x, t) = e^°^(f){x)e-"°^ . (B3) 

In this representation the interaction Hamiltonian reads: 

Hi{t) = e^'>^Hie-^°^ , (B4) 

while 
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Ho{t) = Ho 



(B5) 



The solution of Eq. ( P2|) is well known: 



iS(r) = Tt- exp 



- Hj{T')dr' 
Jo 



Let us evaluate the thermodynamic potential Q: 

From Eqs. and (^7^ we get 

n = -^In Tr (e-^^«5(/3) 



Defining 



-j- In Tr(e-^^°; 



we get 



Whence 



P Tr e-/3^o 



Using Eq. (p6| ) we rewrite (pil| ) as 



One can show that 351: 



/3 

conn ■ 



(B6) 



(B7) 



[B8) 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 



P ^1 ml Jo 
This last equation has been used in Section 8. 

We would like, now, to discuss the thermal propagator of a free scalar field with mass 
m. From the well known expansion: 
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(t>{x, 0) = ^ Y.l'^+'P-^ + e-'P-Vl 



we readily obtain: 



where Ep = ^/p^ + m'^, and we used: 



(B14) 



(B15) 



(B16) 



(B17) 



We are interested in the thermal propagator: 



(r.0(f,r) 0(0) )(' = Gp{x,r) 



(B18) 



Using Eq. (pi5|) and 



P2 



(B19) 



we obtain 



/at a- = —2^ 



Gaix, t) = — 



^ip-x — EpT ^ — ip-x—EpT 

+ 



Now we observe that 



dz 



(B20) 



(B21) 



(B22) 



1 - e-/3^p 7r 27ri (1 - e"^^) (z - ^^p) 

where the integral in the complex z-plane is on the contour F shown in Fig. 16. The 
integrand in Eq. (|B22|) goes to zero exponential when \z\ oo. Thus we deform the 



(B23) 



contour F in F' (see Fig. 16). Applying the Cauchy's integral theorem we get: 



1 - e-^^r> Ep- lUn 
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where ujn = 27mf3. The other term in Eq. (P21|) can be dealt with in a similar way. Thus 
we obtain: 



^ip-X — lUJnT 



-ip-X + iLU„T 



+ 



(B24) 



Finally, observing that y J2p ~^ I performing the change of variables p —p, 

uJn — > ~i^n in the second term in the right hand of Eq. (P24]) , we get: 



d^p e 



■iuj„T+ip-x 



In the following we need to evaluate the sum: 



-E 

P n < + El 

To do this, we use the Sommerfeld- Watson transform [Q: 

dzf{z) + 

1 c+ico+e 1 

^/ [f{z) + f{-z)]-^^A, + A 

Z7^^ J-ioo+e eF'^ — i 



(B25) 



(B26) 



+00 

P n=—oo 



1 

27a 



In our case 



Let us consider, firstly, Ai. We have 



A, 



27ri 



E^ 



^ 

■ioQ Z — 



(B27) 



(B28) 



(B29) 



If r > we close the integration contour in the semiplane Re z < 0, while for r < the 
contour is closed in the semiplane Rez > 0. In this way, by applying the residue theorem 
we obtain: 

e-Ep\r\ 



Ai 



2E„ 



In the same way we get 



A, 



2Ep ei^^p - 1 



Combining Eqs. (|B30| ) and (|B31|) we obtain the desired result: 



P^ul + El 2Ep 



+ 



2E„ 



ei3Ep _ I 



(B30) 



(B31) 



(B32) 
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FIG. 1 Second order corrections to the Gaussian effective potential. 
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FIG. 4 The two-loop (full lines) and the second order (dashed lines) generalized Gaussian effec- 
tive potential for = 1 and three different values of A. 
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FIG. 6 Mass counterterm contributions to the generalized Gaussian effective potential in the 
second order approximation. 
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FIG. 7 The generalized Gaussian effective potential for u = 2 and tfiree different values of A 
with the two-loop corrections. 
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FIG. 8 The generalized Gaussian effective potential for v 
with the three-loop corrections. 



2 and three different values of A 
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FIG. 10 Lowest order thermal correction to the generahzed Gaussian effective potential for v = 1 
and A = 4. Dashed lines refers to the high-temperature expansion. 
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FIG. 12 The critical temperature versus the couphng A for the one-loop effective potential 
(dotted line), the Gaussian effective potential (dashed line), and the generalized Gaussian effective 
potential (full line) in two spatial dimensions. 



61 




FIG. 13 Thermal Feynman diagrams contributing to the second order thermal corrections to 
the generalied Gaussian effective potential. 
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FIG. 14 The second order thermal corrections to the generalized Gaussian effective potential 
for = 1, A = 4 and /3 = 1. 
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FIG. 16 The contours V and V in the complex z-plane. 
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